image 2 of any point on the plane in 
image 1. Any feature pair that is incon- 
sistent with the homography is thrown 
out. The output of the process is a set of 
feature pairs, and the homography. 

The algorithms in this innovation are 
well known, but the new implementa- 
tion improves the process in several 
ways. It runs in real-time at 2 Hz on 64- 
megapixel imagery. The new Shi-Tomasi 
corner detector tries to produce the re- 


quested number of features by automat- 
ically adjusting the minimum distance 
between found features. The homogra- 
phy-finding code now uses an imple- 
mentation of the RAN SAC algorithm 
that adjusts the number of iterations au- 
tomatically to achieve a pre-set probabil- 
ity of missing a set of inliers. The new in- 
terface allows the caller to pass in a set of 
predetermined points in one of the im- 
ages. This allows the ability to track the 


same set of pointsthrough multipleframes. 

This work was done by Daniel S. Clouse, 
Yang Cheng, Adnan I. Ansar, David C. Trotz, 
andCurtisW. Padgett of Caltech for NASA's Jet 
Propulsion Laboratory. Further information is 
contained in a TSP (see page 1). 

The software used in this innovation is 
available for commercial licensing. Please con- 
tact Daniel Broderick of theCalifornia Institute 
of Technology at danidb@caltech.edu. Refer to 
NPO-46916. 


@ Sparse Superpixel U nmixing for H yperspectral I mage Analysis 

NASA's Jet Propulsion Laboratory, Pasadena, California 



Here a Subwindow of Observation demonstrates superpixel segmentation. Left: original subimage. 
Center: coarse segmentation, minimum region size 100. Right: fine segmentation, minimum region size 20. 


Software was developed that automat- 
ically detects minerals that are present 
in each pixel of a hyperspectral image. 
An algorithm based on sparse spectral 
unmixing with Bayesian Positive Source 
Separation is used to produce mineral 
abundance mapsfrom hyperspectral im- 
ages. A "superpixel" segmentation strat- 
egy enables efficient unmixing in an in- 
teractive session. 

The algorithm computes statistically 
likely combinations of constituents 
based on a set of possible constituent 
minerals whose abundances are uncer- 
tain. A library of source spectra from 
laboratory experiments or previous re- 
mote observations is used. A superpixel 
segmentation strategy improves analy- 
sis time by orders of magnitude, per- 
mitting incorporation into an interac- 
tive user session (see figure). 

Mineralogical search strategies can be 
categorized as "supervised" or "unsuper- 
vised." Supervised methods use a detec- 


tion function, developed on previous 
data by hand or statistical techniques, to 
identify one or more specific target sig- 
nals. Purely unsupervised results are not 
always physically meaningful, and may ig- 
nore subtle or localized mineralogy since 
they aim to minimize reconstruction 
error over the entire image. This algo- 
rithm offers advantages of both methods, 
providing meaningful physical interpre- 
tations and sensitivity to subtle or unex- 


pected minerals. 

This work was done by Rebecca Castano 
and David R. Thompson of Caltech and 
Martha Gilmore of Wesleyan University for 
NASA's Jet Propulsion Laboratory. For more 
information, contact iaoffice@jpl.nasa.gov. 

The software used in this innovation is 
availablefor commercial licensing. Pleasecon- 
tact Danid Broderick of theCalifornia Insti- 
tute of Technology at danielb@caltech.edu. 
Refer to NPO-47038. 


@ Intelligent Patching of Conceptual Geometry for CFD Analysis 

Langley Research Center, H ampton, Virginia 


The iPatch computer code for intelli- 
gently patching surface grids was devel- 
oped to convert conceptual geometry to 
computational fluid dynamics (CFD) 
geometry (see figure). It automatically 
uses bicubic B-splines to extrapolate (if 
necessary) each surface in a conceptual 
geometry so that all the independently 
defined geometric components (such as 
wing and fuselage) can be intersected to 
form a watertight CFD geometry. The 
software also computes the intersection 
curves of surface patches at any resolu- 
tion (up to 10 -4 accuracy) specified by 


the user, and it writes the B-spline sur- 
face patches, and the corresponding 
boundary points, for the watertight CFD 
geometry in the format that can be di- 
rectly used by the grid generation tool 
VGRID. 

iPatch requires that input geometry 
be in PLOT3D format where each com- 
ponent surface is defined by a rectangu- 
lar grid -R x( i ,j ) , y(i,j), z( i,j) ) : 1 < i < m, 1 
< j < n}that represents a smooth B-spline 
surface. All surfaces in the PLOT3D file 
conceptually represent a watertight 
geometry of components of an aircraft 


on the half-space y> 0. Overlapping sur- 
faces are not allowed, but could be fixed 
by a utility code "fixp3d". The fixp3d 
utility code first finds the two grid lines 
on the two surface grids that are closest 
to each other in Hausdorff distance (a 
metric to measure the discrepancies of 
two sets); then uses one of the grid lines 
as the transition line, extending grid 
lines on one grid to the other grid to 
form a merged grid. 

Any two connecting surfaces shall 
have a "visually" common boundary 
curve, or can be described by an inter- 
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